Witnessing Macroscopic Entanglement in a Staggered Magnetic Field 
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We investigate macroscopic entanglement in an infinite XX spin-^ chain with staggered magnetic 
field, Bi — B + e~' lnl b. Using single-site entropy and by constructing an entanglement witness, 
we search for the existence of entanglement when the system is at absolute zero, as well as in 
thermal equilibrium. Although the role of the alternating magnetic field b is, in general, to suppress 
entanglement as do B and T, we find that when T = 0, introducing b allows the existence of 
entanglement even when the uniform magnetic field B is arbitrarily large. We find that the region 
and the amount of entanglement in the spin chain can be enhanced by a staggered magnetic field. 
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Quantum entanglement is a fundamental aspect of 
quantum physics. It demonstrates the non-local nature 
of the theory in that an entangled system contains corre- 
lations that cannot be described by its subsystems alone. 
Instead these quantum correlations are attributed to the 
overall system Further, entanglement is an impor- 
tant resource in quantum information and computation. 
In particular, solid state quantum computation has be- 
come a topic of much research and several proposals for 
physical implementation have been investigated. The 
Heisenberg interaction is the model used in many phys- 
ical applications of quantum computation, for example, 
quantum dots [2] and cavity QED It has also been 
shown that the Heisenberg interaction can be used to im- 
plement any circuit required by a quantum computer Q • 
Therefore, entanglement in one-dimensional spin chains 
has been the subject of much interest. This entangle- 
ment has been studied both in the case of a finite spin 
chain 5, 6] and in the thermodynamic limit 7] where the 
length of the spin chain becomes infinite. 

Macroscopic entanglement is a more recent concept. It 
demonstrates that non-local correlations persist even in 
the thermodynamic limit. This type of entanglement can 
be detected by measuring macroscopic quantities such 
as internal energy and magnetic susceptibility [l3| as it 
has been proven that such quantities can be used as en- 
tanglement witnesses. It has been shown experimentally 
[13, [TH, that the behaviour of observable macroscopic 
quantities such as magnetic susceptibility depends, most 
significantly at low temperatures, on entanglement. This 
demonstrates that entanglement is vital in the explana- 
tion of how macroscopic materials behave. Macroscopic 
entanglement in a Heisenberg spin chain has been stud- 
ied previously Q only for a uniform magnetic field. The 
Hamiltonian of this chain is used to construct an entan- 
glement witness 8, 9, 10] which shows that entanglement 
disappears for high uniform magnetic field just as it does 
for high temperature. 

In real systems, the magnetic field need not be the 
same at each site in the chain. In solid state systems, 
there exists a possibility that an inhomogeneous Zee- 



man coupling could induce a non-uniform magnetic field. 
Moreover, an experimental system is likel y t o contain 
magnetic impurities. Copper Benzoate [Tl|] , [12j is a prac- 
tical example of a system in a non-uniform magnetic field. 
In this case, the alternating field is in a direction per- 
pendicular to the uniform field. Alternatively, such im- 
purities could be introduced artificially. Therefore, the 
possibility that such a field could affect entanglement, 
whether to reduce or increase it, is an important subject 
to investigate. In reality, systems have a finite temper- 
ature so the thermal case must be considered. Hence, 
in this paper, we discuss the effect of a site dependent 
magnetic field on thermal macroscopic entanglement in 
a 1-D infinite spin-| chain. We also consider the zero 
temperature case. Interestingly, we show that an alter- 
nating magnetic field can compensate for the effect of a 
uniform magnetic field at T = 0. 
The Hamiltonian considered is 
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where J is the coupling strength between sites, and 



B, = B + e 



— ittI 



b is the site dependent magnetic field. 



Although such a field is not likely to occur in nature, our 
work allows us to investigate how a non-uniform magnetic 
field affects entanglement in a model which is analytically 
solvable. A similar Hamiltonian with cyclic boundary 
conditions has previously been diagonalized [l6| using a 
method first set out by S. Katsura [13, EH- In our dis- 
cussions of a finite JV-spin chain, we consider the case of 
open boundary conditions with N even. In fact, these 
constraints are no longer relevant in the thermodynamic 
limit and hence our conclusion is unchanged if, for exam- 
ple, the cyclic model is used. 

To identify entanglement in this system, we use an 
entanglement witness, i.e. an operator whose expectation 
value is bounded for any separable state. The power of 
our witness is such that we can identify the existence of 
entanglement even for a thermal system which is a mixed 
state in general. Alternatively, single-site entropy can be 
taken as evidence of entanglement when T = since 



the total system is in a pure state. The purity of the 
single-site density matrix shows that the entanglement 
witness is not optimal at absolute zero. Thus we find 
that although the alternating magnetic field, b, acts in 
general to suppress entanglement similarly to B and T, 
at zero temperature, increasing b allows the system to be 
entangled for arbitrarily large B. Hence the effect of the 
staggered field at T — is to increase both the amount 
and the region of entanglement. 

Entanglement - A system is said to be entangled when 
its density matrix cannot be written as a convex sum 
of product states. For a pure state, dividing the system 
into two subsystems A and B allows the von Neumann 
entropy to be used as a measure of entanglement. If we 
trace section B out of the density matrix to find pa, the 
von Neumann entropy, S(pa) — — Tr(p^log 2 pa), can be 
calculated. S(pa) — corresponds to a separable state 
while when S(pa) — 1, the system is maximally entan- 
gled. In the case of a mixed state, there is no unique mea- 
sure of entanglement for a multipartite system. However, 
we can construct an entanglement witness. 

An entanglement witness is an operator whose ex- 
pectation value for any separable state is bounded by a 
value corresponding to a hyperplane in the space of den- 
sity matrices. An entanglement witness is only a suffi- 
cient condition for the existence of entanglement. Hence 
failure of the witness to detect entanglement does not 
necessarily mean the system is separable. Though wit- 
nesses simply detect rather than give a measure of en- 
tanglement, they have significant advantages over other 
methods. For example, they naturally incorporate tem- 
perature and many witnesses, such as magnetic suscep- 
tibility, can be experimentally measured [la ]. 

The partition function and entanglement witness - 
Many thermodynamic variables can be derived from the 
Helmholtz free energy, F = — TlnZ, where Z is the par- 
tition function. As dF/dX — we see that when 
X = B, we obtain the magnetization M — J2i( a i) — 
—dF/dB. In particular, we define the entanglement wit- 
ness 



W = 



2 d\nZ 1 



PN dJ 



N ^ 
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(2) 



where (3 — 1/T. Our witness, W, identifies a larger en- 
tangled region than witnesses found previously ,7]. In a 
separable state, it satisfies the bound \ W\ < 1, which can 
be shown as follows. With p = J2i Pi Pi ®P2®'"®Pn^ we 
have \{afaf +1 ) + {afaf +1 )\ = | (^f)K +1 ) + K)K +1 )| < 



T y\2 



'l I -r \"i I y \"f+i> 2 + W+i) 2 < 1 for an y l - Thc 
upper bound for the inequality is found by using the 
Cauchy-Schwarz inequality and the condition that for 
any state, (erf) 2 + (of) 2 + (°f ) 2 < 1- Thus, any state 
that violates the inequality \W\ < 1 is entangled. 

Diagonalization of Hamiltonian - In order to find the 
partition function of the system, we must diagonalizc 
the Hamiltonian, Eq. (TT]) and find its energy eigenval- 
ues. The open ended Hamiltonian can be exactly diag- 




FIG. 1: The surface \W\ = 1 in the space of magnetic field, 
B, temperature, T, and alternating magnetic field, b. The 
region between the surface and the axes is entangled. 



onalized in a standard way via several steps: a Jordan- 
Wigner tansformation, a Fourier transformation and fi- 
nally a Bogoliubov transformation. The Jordan- Wigner 
transformation, 



at 



i-i 

n< 

m — 1 



ipf 



ia. 



(3) 



maps the Pauli spin operators into fermionic annihilation 
and creation operators ai and a\. These satisfy the anti- 
commutation relations {et/,afc} = and {a/,a[,} = 5;^. 
Preserving the anti-commutation relations, the operators 
can now be transformed unitarily using a Fourier trans- 
formation, ai = J2k=i d k Bin(-j^) and by a Bo- 
goliubov transformation, 

dk = ak cos 9k + (3k sin 9k 
dN+i-k = Pk cos 9 k - a k sin 9 k . 



(4) 



Setting tan26>fc = b / [J cos(ir k/(N + 1))] eliminates the 
off-diagonal terms leaving the Hamiltonian in diagonal 
form 



H = 



N/2 

E 

k=l 



(5) 



where = 2B ± 2^J 2 cos 2 (7rfc/(A + 1)) + b 2 . Thc 
operators oth and (3k satisfy the anti-commutation re- 
lations {a k ,a\} = {(3 k , Pj} = Skd and {a k ,(3i} = 
{ak,(3j} — 0. Using the eigenvalues of the Hamilto- 
nian, we find that the partition function can be written 

z = Uk=i 2 cosh (/ 3A fc I 2 ) 2 cosh (f j K I 2 ) ■ In the ther- 
modynamic limit, N — > oo, we can treat ui = irk/N as a 

continuous variable, to find 



N 

InZ = — 

7T 



tt/2 



duj In 



4 cosh 



cosh 



(PK 



\ 2 



(6) 



3 



where = 2B ± 2 V 'J 2 cos 2 Hi + b 2 . We can now use Eq. 
@ to calculate the entanglement witness for our system. 

The region of entanglement detected by our witness 
has been plotted in Fig. JT]). The figure shows the re- 
gion of uniform magnetic field, B, temperature, T, and 
alternating magnetic field, b, within which we always find 
entanglement. At fixed values of b, the entangled region 
of the T — B plane shrinks as b increases until a critical 
value is reached above which entanglement is no longer 
detected. Consider now a plane perpendicular to the T 
axis. At zero temperature we see that until the critical 
value b c = 0.56, increasing b has no effect on the value 
that the uniform magnetic field can take with the system 
remaining in an entangled state. Above b c , our witness 
detects no entanglement. However, we later show that 
this witness is not optimal at zero temperature. 

Our witness shows entanglement behaving as we would 
expect physically. A high temperature causes the system 
to become mixed. Thus no quantum correlations can 
survive and the system is separable. Moreover, the effect 
of the magnetic fields is understandable as local oper- 
ations which enhance classical correlations. When the 
uniform magnetic field becomes large, spins tend to line 
up in a direction parallel to that field. This clearly de- 
creases quantum correlations in the system. Using the 
same reasoning, if the alternating magnetic field is large, 
the spins tend to anti-align which is also a product state. 
Hence, as Fig. ([IJ shows, all of these parameters cause 
the system to become separable if they are large enough. 
Interestingly, we find a counter example of this expected 
behaviour and identify a region where the system is en- 
tangled even in a large magnetic field. We discuss this in 
the following section. 

Single-site entropy - When T = 0, the total system is 
in a pure state, i.e. the ground state. Thus, if the density 
matrix of a single spin is in a mixed state, the particle 
must be entangled with the rest of the spins. This can 
be quantified using the entropy of a single spin. The 
Z-th spin density matrix, pi = | J2ie{i, x ,y,z} cr l( cr l)' can 
be obtained from the total density matrix p = e~° H /Z. 
Moreover, we find that the single-site density matrix is 
readily diagonalized since (erf) = (af) = 0. This follows 
from the fact that erf and af are linear combinations of 
the fermion operators a>k, oi\ 1 (3k, and f3\, all of which 
have zero expectation values. 

In the thermodynamic limit N — > oo, the system is 
translationally invariant for all odd sites and for all even 
sites. Hence the single-site magnetization (erf) can be 
obtained from the total magnetization M and the total 
staggered magnetization M s . In the limit of zero temper- 
ature, these are given by 

v^, „ 1 d ( 2fl\ 

M = Y J {«l) = ~^Z = N(l--) (7) 

1 f) AT ft* 



Even Site Odd Site 




FIG. 2: We plot the entropy for an even site spin and an odd 
site spin respectively when T — and J = 1. Excluding the 
region B > \j 1 + b 2 , the single-site entropy is non zero and 
each spin is entangled with the rest of the system. The odd 
site entropy shows a maximum when B = b + e. 

where f(uj) — 26/ V 'J 2 cos 2 to + b 2 and 

!7r/2 for B < b 

cos- 1 (^/B 2 -b 2 /J^j for b < B < v 7 ./ 2 + b 2 
for B > V J 2 + b 2 

To obtain these results, we have used 
lim^^oo tanh(/3a;) = af/|af|. By virtue of translational 
invariance, (erf), is the same for all even sites and for all 
odd sites. Thus defining (erf) = Tr(erf p) = 2r\i — I, such 
that the single-site density matrix is pi — diag(?7;, 1 — 77;), 
we have 

m = l{l + jf[M+(-l) l M B ]}. (8) 

From this, we can obtain the entropy of the Z-tli spin, 
S = -nt log 2 rji - (1 - 77;) log 2 (l - 77/). 

The single-site entropy for odd and even sites is plotted 
in Fig. @ with J = 1. In both cases, we find entangle- 
ment for various values of B and b even when the mag- 
netic fields are large. The entropy, and therefore entan- 
glement between each spin and the remainder of chain, is 
non-zero everywhere except when B > V 'J 2 + b 2 . Hence 
entanglement exists when the coupling strength between 
nearest neighbour spins, J, is more than \/B 2 — b 2 . We 
note that this corresponds to when the square of the in- 
teraction strength is greater than the product of the to- 
tal magnetic field on two adjacent sites. In addition, 
we observe from Fig. ^ that the maximum single- 
site entropy occurs when both magnetic fields B and b 
are zero. Introducing the magnetic fields reduces the 
amount of entanglement in the system except along the 
peak in the odd site entropy. We find from the max- 
imum entropy, Si = 1, that the maximum entangle- 
ment occurs in the region b < B < V J 2 + b 2 when 

- 7r/2 = (-1)' J n bj\J J 2 cos 2 to + b 2 duj is satisfied. 
This corresponds to when the magnetization, M, is equal 
to the staggered magnetization, M s . For even sites, there 
is only one solution to this at B = b = 0. For odd sites 
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FIG. 3: Plot of the maximum entropy Si = 1 for odd site spins 
with different values of J. As the value of J increases, this 
maximum entropy moves further away from the line B — b. 

however, this is satisfied for any finite uniform magnetic 
field B at B = b + e where £ is a positive value. The so- 
lutions in Fig. ([3]) correspond to the peak in the odd site 
entropy and so occur within the range b<B< VJ 2 + b 2 . 
In general, s becomes larger as J increases as shown in 
Fig. ([3]), and becomes smaller as B and b increase. We 
note that the peak does not occur at B = b. As B 
and/or b tend to infinity, the amount of entanglement 
in the system tends to zero. Further, as the magnetic 
fields increase, the curves in Fig. ([3]) tend to the B = b 
line. At B = b, the system is no longer maximally en- 
tangled. At this point, we find from the Hamiltonian 
that odd sites have zero magnetic field, and even sites 
have field strength 2B. Hence at the peak, odd site spins 
have a small magnetic field, e, while in comparison, even 
site spins have a large field, 2B — e. The implications of 
this peak are that even in the limit of large (though not 
infinite) uniform field £?, an odd site can be maximally 
entangled with the rest of the system if an appropriate 
alternating magnetic field, b is introduced. In practice, 
fine-tuning b to achieve maximal entanglement may be 
difficult. However, for any £?, sufficiently increasing b 
(when b > \f B 2 — J 2 ) will create entanglement in the 
chain. Hence a staggered magnetic field can enhance the 
amount of entanglement present in the system. 

Although the single-site entanglement relates only to 



the zero temperature case, changing the temperature by 
a small amount should not change its behaviour. Hence 
for very low temperatures, we see that the entanglement 
witness is not optimal. 

Conclusions - Our best estimate for finite temperature 
entanglement is the witness which shows b reduces the 
region of entanglement in the chain. That is, a staggered 
magnetic field reduces the entangled region. If this be- 
haviour is true even for an optimal finite temperature 
witness, these results have consequences for larger scale 
quantum computation in solid state systems. As inhomo- 
geneities in the magnetic field exist naturally in the Zcc- 
man coupling between atoms, the region of entanglement 
is naturally decreased compared to when 6 = 0. Quan- 
tum computation relies on entanglement so as introduc- 
ing b reduces both the temperature and uniform magnetic 
field at which entanglement persists, our result shows it 
may be more difficult than previously thought to con- 
struct useful quantum computers using one-dimensional 
systems. 

Our entanglement witness is invaluable as by applying 
it to our system, it allows us to see how temperature af- 
fects the entanglement in the spin chain. However, the 
witness does not tell us how the entanglement actually 
behaves in the presence of B, T or b as it does not detect 
all entanglement in the system. Conversely, the single- 
site entropy shows us exactly how the entanglement is 
affected by the magnetic fields at zero temperature, al- 
though it is unknown how to extend this entropy to a 
finite temperature. Using this entropy, we have shown 
that in the thermodynamic limit, a staggered magnetic 
field enhances both the region and the amount of entan- 
glement in our spin chain. Hence, both the witness and 
the entropy are essential in characterizing the entangle- 
ment in the system. Further, the region of entanglement 
identified by the witness is consistent with that of the 
single-site entropy. If the behaviour of the entanglement 
as shown by the entropy persists at higher temperatures, 
we may be able to counteract any Zeeman coupling by ap- 
plying an appropriate magnetic field, hence maximizing 
entanglement for odd sites. This will be an interesting 
topic for future research. 
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